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ABSTRACT 

We present a new method of studying quadruple lenses in elliptical power-law 
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^ . potentials parameterized by ^(x,y) oc (a;^ + y'^ Iq^Y/"^ I (3 (0 < /? < 2). For 

CO ' this potential, the moments of the four image positions weighted by signed 

f^ ' magnifications (magnification times parity) have very simple properties. In 

<3\'. 

Qs , particular, we find that the zeroth moment - the sum of four signed magnifi- 

(~| . cations satisfies ~ 2/(2-/3); the relation is exact for /? = (point-lens) and 

I ' /3 = 1 (isothermal potential), independent of the axial ratio. Similar relations 

o ■ 

Vh ' can be derived when a shear is present along the major or minor axes. These 

jrt ! relations, however, do not hold well for the closely-related elliptical density 

distributions. For a singular isothermal elliptical density distribution without 
j^ ' shear, the sum of signed magnifications for quadruple lenses is w 2.8, again 

C^ ' nearly independent of the ellipticity. For the same distribution with shear, 

the total signed magnification is around 2-3 for most cases, but can be sig- 
nificantly different for some combinations of the axial ratio and shear where 
more than four images can appear. 
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1 INTRODUCTION 

Since the discovery of the first double lens 0957+0561 by Walsh, Carswell & Weymann 
(1979), gravitational lensing has found many cosmological applications (see Schneider, Ehlers 
& Falco 1992; Narayan & Bartelmann 1998 and references therein). One of these applica- 
tions is to probe the potential of lensing galaxies in multiply-imaged quasars. In this aspect, 
the quadruple lenses are ideal systems since they provide more constraints than the dou- 
ble lenses. Previous modelings of these systems use mainly isothermal elliptical potentials 
and their variants (e.g., Kochanek 1991; Keeton, Kochanek & Seljak 1997; Witt, Mao & 
Schechter 1995). In many of these modelings, attention was paid only to the image posi- 
tions with a few exceptions, which considered the fiux ratios (e.g., Keeton, Kochanek & 
Seljak 1997). Since optical fiux ratios are likely to be infiuenced by microlensing (Chang 
& Refsdal 1979), they should be used with care. On the other hand, radio fiuxes should 
be nearly unaffected by microlensing since the radio emission regions are believed to have 
larger sizes (there are exceptions for very highly magnified systems such as B1422+231, 
see Mao & Schneider 1998). Satisfactory lens models should ideally fit the image positions 
as well as the fiux ratios. However, previous numerical studies often find it very difficult 
to fit these with isothermal potentials and their variants (e.g., Keeton, Kochanek & Seljak 
1997). The purpose of this paper is two- fold: first we extend the isothermal potentials to a 
more general class of power-law potentials and uncover some generic scaling relations. The 
three-dimensional correspondence of this class of potential provide a reasonable analytical 
approximation to the galactic potential (Evans 1994). Second, we study the relations for 
related isothermal elliptical density distributions and show that there are qualitative and 
quantitative differences between elliptical density and elliptical potential distributions. 

The outline of the paper is as follows. In §2, we introduce the elliptical power-law po- 
tentials and derive some scaling relations with a new moment approach. In §3, we study 
the singular isothermal elliptical density distribution with arbitrary shear to understand 
the difference between these two types of potentials. In §4, we summarize and discuss the 
implications of our results for lens modeling. 

2 ELLIPTICAL POWER-LAW POTENTIAL 

In this section we study elliptical power-law potentials as models for lensing galaxies. This 
potential form is widely used in gravitational lensing (e.g., Blandford & Kochanek 1987; 
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Kassiola & Kovner 1993 and references therein). It is easier to study than the more reahstic 
but more complex elhptical density distribution (e.g., Kormann, Schneider & Bartelmann 
1994; Hogg & Blandford 1994; Chae, Khersonsky, & Turnshek 1998); we return to the 
elhptical density distributions in §3. 
The power-law potential is given by 

^{x,y) = a ^ ^y^ ^ for < /3 < 2 (1) 

where g (0 < g < 1) is the axial ratio of the lensing galaxy, (3 the slope of the potential 
and a determines the overall strength of the lens system. The coordinate system has been 
chosen such that the center of the lensing galaxy is at the origin and the major axis of 
the galaxy is along the x-axis. In general the orientation of the major axis of the lensing 
galaxy is unknown, but can be easily determined if the image configurations can be modeled 
by an elliptical potential (Witt 1996). Note that the point-mass potential is described by 
/? = 0, g = 1, the elliptical isothermal sphere by /? = 1, and the constant surface mass 
distribution by /3 -^ 2, g = 1. Thus the commonly used potentials in numerical modeling are 
all special cases of this class. 

2.1 Lens equation 

The lens equation for the power-law potential is given by 

^ = ^~ (x2+yVg2) 1-/3/2' (^) 

^yr' (4) 



^ (X2 + 2/2/^2)1-/3/2- 

It is easy to show that the lens equation with an on-axis shear can be transformed into 
eqs. (^ and @) with suitable coordinate transformations as given in §2.5. So with small 
modifications our results apply to the case with an on-axis shear as well (§2.5). Notice that 
the elliptical potentials (including the special power-law form studied here) restrict the image 
and galaxy positions to be located on a hyperbolic line (Witt 1996). Image configurations 
which do not follow these constraints require far more complicated lens models, in particular 
may require an off- axis shear (Keeton, Kochanek & Seljak 1997; Witt & Mao 1997). 
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2.2 Exact relations for /? = 1 and P = 

Although the lens equations (|^) and (Q) are in general too complicated to solve for the image 
positions and the magnifications, here we present a new method to directly relate the image 
magnification and positions to the properties of the galaxy potential. The new method is 
similar to the Jeans equation in galactic dynamics (e.g., Binney & Tremaine 1987, p. 195). 
It is possible to extend this new technique to other potential forms such as that for binary 
lens. However, in the latter case, the application may be limited since in microlensing we 
can not resolve the individual images. In the next subsection, we illustrate our technique 
using an isothermal {j3 = 1) elliptical potential; in §2.3, we discuss the case for /5 = 0. 

2.2.1 Isothermal elliptical potential ((3 = 1) 

There are different ways to derive relationships between the image positions and their mag- 
nifications. Sometimes it is more convenient to use the complex formalism (cf. Rhie 1997 for 
the binary lens application) or the resultant method (e.g., Witt & Mao 1995; Dalai 1998). 
Here we adopt a general analytical approach using real quantities. When a source is inside 
the inner caustic of an elliptical potential, there are always four images. The lens equation 
(0) or (^) can be manipulated into two polynomials of degree 4 of separated arguments in x 
and y: 

4 4 4 4 

9x{x,^,r]) = J2(^nx'' = a4Y[{x-Xi), gy{y,^,v) = Y.^^y'' = ^^IKv ~ Vi) (5) 

n=0 «=1 n=0 j=l 

where (xj, yi) are the coordinates of the i-th image, and the polynomial coefficients are given 
by 

a, = (1-gY, b, = q\l-q'f, (6) 

as = 2e(l-g^)(g^-2), 63 = 2w^(l - g^)(2g2 - 1), (7) 

a, = -a\l-qy + 6e{l-q') + q\q'e + v') h = a^ - 6{vY + e){l - q') (8) 

a, = 2^[a'{l - q') - 2^^ + q'{e - v')[ h = -2qMa'il - q') + 2gV + q'ie - v')] (9) 
ao = -eia' - ie + qV)], bo = -vY[a' - (e + q'v')] (10) 

The Jacobian matrix for the mapping from the lens plane to the image plane is defined as 

^ dx dy 

For each image, there is a one-to-one inverse mapping from the image plane to the lens 
plane, 
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( dx dx 

d^ dri 

dy dy 

d^ dri 



-*- / dy dy 

^ dx dx 



(12) 



where the second equahty comes from the conventional definition of an inverse matrix. A 
comparison of the diagonal elements gives 
dy I d^ dx 1 Or] 



dri det J dx' d^ det J dy 
The sum of the two expressions is 



dy dx 



d^ dri 
dx dy 



drj d^ det J 
where we have used 



|^ + ^ = l + detJ for 13 = 1. 
dx dy 



det J 



+ 1, 



(13) 



(14) 



(15) 



(16) 



The signed magnification for a given image is by definition 
1 
det J (x„j/,) 
where /Xj is the absolute magnification for a given image, pj = ±1 is the parity of of the 

images, and we substituted det J from eq. ([T^). 



% dxj 
drj d^ 



Now the sum of signed magnifications is then 



d 



d 



d 



ag 



d 



(17) 



In the second step we have used the fact that the sum of the roots of a polynomial is related 



(18) 



to its coefficients. Using the coefficients as given in eqs. (|]) to (pJOf ), we obtain 

4 

i=l 

A similar relation was found recently by Dalai (1998) for (3 = 1. This implies that for an 
isothermal elliptical potential the signed magnification is always 2, independent of the axial 
ratio and the image configurations. 

Using similar techniques, we can calculate higher order moments of positions, defined as 



Mn,x = "YpifiiXl', Mn,y = ^Pi^-iyl 



i=l 



i=l 



For illustration, we derive the x-component of the first moment: 



■^l,x = J2Pil^i^i = E 



j=l 



j=l 



dyi dxj 
drj d^ 



Xi 



(19) 



(20) 



where in the second equality we have again used eq. ([1^) . Changing now the order of sum- 
mation and derivation, and using the fact that the Jacobian matrix is symmetric {dy/dC, = 
dx/drj), we obtain: 
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4 I d ^ d '^ ^ 



(21) 



For the first and third terms we can easily replace the sum by expressions of the coefficients 
of the polynomial: 



E 



X, 



2a2a4 



Y.yl 



hi - 26,6, 



'2(^4 



/)2 
"4 



and E-^« ~ • 

i=l «4 



(22) 



i=l "4 i=l 

For the second mixed term in x and y of eq. (0), we recall that for any elliptical potential, 
the X and y components for a given image are related (cf. eqs. and 0]) by: 



rjq^Xi 



= g — , thus yi = - 

V-Vi Vi ^-Xi{l 

Now the mixed term can be written as 



q' 



(23) 



^XiVi 



rjq 



E 



X7 



i=l 



^-Q'tl ^/(l - '?') - ^^ 



rjq 



l-g2 



4e 



e^ 



E 



1 



(24) 



03 

where we have expressed the mixed term in terms of a sum of different powers of Xj. The 
last term in the previous equation can be expressed in terms of the polynomial: 

' 1 , 1 dg,{x,i,Ti) 2(l-g2^ 



E 






=«/(i"«^) g,{x,i,rj) 



dx 



iq' 



Collecting all terms together, we obtain 

4 4 

^Pi/iiX, = 2^, Yp,iJ,iyi = 2r], 



(25) 



(26) 



i=l 



i=l 



a remarkably simple result. Higher order moments can be obtained similarly and are sum- 
marized in §2.4. However, first we discuss results for another special case with /? = 0. 



2.3 Point lens {(3 = 0) 

For (3 = 0, the potential approximates a point lens plus shear. It is easy to separate the lens 
equations into two polynomials of degree 4 for x and y - We only have to insert eq. (p3D 
into the lens equations (^) and (1). Therefore we do not present the coefficients here. The 



derivations for magnification moments are similar, except that eq. (|1^) must be changed. To 
do this, we notice that for /5 = 0, we have 



dx dy 



1 + g'. 



Using eq. (p!3D, we obtain 



dy 



,dx 



1 



-^ + q — 

dr] d^ det J 



dx dy 



1 + g^ 
det J 



(27) 
(28) 
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In the last step, we have expressed the signed magnification {fiiPi = 1/det J) in terms of 
dy/drj and dx/dC,, similar to the case for (3 = 1. The rest of derivation proceeds nearly 
identically. For the zeroth and first moments we have 

4 4 4 

^/iiPi = i, ^piiiiXi = i, Y.Pif^iyi = v, (29) 

i=l i=l i=l 

2.4 Relations for more general power-law potentials 

For more general cases with non-integer values between and 2, there is no simple way 
to derive analogous relations. However, the special cases for /3 = and (3 = 1 suggests 
approximation by interpolation. In particular, 

X:/i.P.^i?, B^j^. (30) 

For /5 = and /5 = 1, eq. (^) is exact. For other non- integer values, this equation is only 
approximate, but we find numerically that it is valid to ^5%. Therefore, for an quadruple 
lens well-modeled by power-law elliptical potentials, this relation is approximately valid, 
independent of the source position. A similar exact relation is also satisfied by the binary 
lens, where the sum of signed magnifications of five image cases is equal to 2 (Witt & Mao 
1995; Rhie 1997). 

Starting from eq. (pO]), we can obtain moment equations that relate the image positions 
and magnifications to the parameters of the potential. 



^liiPiXi = Bi, ^^iiPiyi = Bri, (31) 

i=l j=l 

4 4 

Y.^i^P^x1 = -A + Bi\ Y.^'^P^y^^=M~^'' + Br,\ (32) 

i=l i=l 

f:/x.p.xf = -Ai{B + \^)+Be, J2f^^P^yf = ^m-'''{B + ^y^) + Brj', (33) 

where A = a^Bq^ /{I — q^). Similar to eq. (pOD, these relations are exact only for /5 = and 
(3 = 1. However, numerical tests show that these formulae are accurate to ^1% for other 
non- integer values of (3 (or B). 

2.5 The inclusion of an on-axis shear 

The lens equation with an on-axis shear is given by 

^ = ^ + ^^- (a:^ + ,Vg2)i-/./2 (34) 
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We can show that these equations can be transformed into equations (^ and (^) by a simple 
transformation. It can be easily verified that the transformation is given by 

r = -4= V'=^== (36) 



x' = xVl + 7 y' = yJl-i (37) 



a = a{l + j)-^/' P' = P q' = q^]^, (38) 

where we denote the transformed quantities with a prime. Applying the above transforma- 
tion, we obtain 

/ _ / a'y'q'-^ 

"^ y (a;/2 + ^/2/g/2)i-/3'/2- y"^^) 

These are identical to eqs. (§) and @ if we drop the prime. 
The magnification transforms like 

K = /^.(l-y) (41) 

Now it is straightforward to generalize the moment relations (eqs. [^ 33]) to the case 
with an on-axis shear by applying eqs. (|36|)- (|4l|) . In particular, for the zeroth and first 
moments we have 

^^^^^^^-(2 -/?)(! -7^) (1 + 7)' t/'^'^' (2-/5)(l-7^)(l-7)' ^ ^ 

Notice that the inclusion of an on-axis shear changes the magnification relation (eq. [^ ) 

only quadratically. It is straightforward but tedious to show that this relation holds even for 

an off-axis shear (cf. Dalai 1998). For the observed gravitational lenses, the external shear 

is usually between < 7 < 0.3 (Keeton, Kochanek & Seljak 1997), this would mean that 

for given slope of the elliptical power-law potential, the total sum of signed magnifications 

is known to within 10%. 

Two interesting observations are worth making here. First, these simple formulae indicate 

that the moments of x (y) components Aix,n of the images are only related to the ^ (77) 

component of the source position. This seems remarkable since the magnification is in general 

a function of both ^ and rj. Second, for an elliptical potential with an on-axis shear, we have 
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7 unknown quantities (^, r], a, q, P, /ii, 7) with 7 moment equations, so in principle we can use 
the moments to solve for the lens parameters directly. Indeed we have implemented such a 
routine. Unfortunately, using Monte Carlo simulations, we find that the results depend rather 
sensitively on the accuracy of image positions and flux ratios. The present observational 
uncertainties do not allow us to carry out such an exercise, although in the future such an 
effort may be more worthwhile in order to determine whether elliptical potentials provide 
adequate fits to the observed lenses. In any case, we will show next that the magnification 
relations do not hold exactly for the elliptical density distributions and therefore these 
relations should be used with some caution. 

3 SINGULAR ISOTHERMAL ELLIPTICAL DENSITY DISTRIBUTIONS 

It is well known that power-law potentials develop unphysical features, such as dumb-bell 
shaped density contours or even negative densities, when the axial ratio is small (e.g., Kassi- 
ola & Kovner 1993). So it is important to check whether the relations derived for the simpler 
elliptical potentials are valid for the more realistic elliptical density distributions. Relations 
concerning only image positions seem to be equally well satisfied by the elliptical potentials 
and elliptical density distributions (Witt & Mao 1997). However, this is not necessarily true 
for the magnification relations since magnifications involve second order derivatives of the 
lensing potential rather than the first order for the image positions, and hence they could 
be more sensitive to the detailed potential forms. We show that indeed elliptical density and 
elliptical potentials have somewhat different magnification relations. In this paper, we focus 
on the singular elliptical density distribution (including an arbitrary shear) for two reasons: 
the first is that this model can reproduce the observed image configurations well, particularly 
the models incorporating a shear (Keeton, Kochanek & Seljak 1997). The second reason is 
that relatively little is known about the magnification properties for the elliptical density 
distributions, so a complete treatment of the general case is beyond the scope of this paper. 
In polar coordinates (r, 9) , the surface density for the singular isothermal distribution is 
given by 

"''•■'> = 1^ (1 -,cL 29) V^ - (^^' 

where e = (1 — q'^)/{l + q"^), q is the axial ratio and a is again a characteristic radius that 
determines the overall image separations. For e > 0, the major axis is along the x-axis. The 
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lens equation (including a shear) is 

^ /97 f-r,c ft 

(45) 



^ = (l + 7i)a; + 72y- ^tan 



ecos2^)V2 
/2esin^ 



V2e [(1 — ecoszb')^/^] 

where (71 = 7 cos 20s, 72 = 7 sin 20s with 7 being the strength of shear and 0s the "direction" 
of shear. The magnification for a given image is 

(det J)-^ = 1-2k{1 + 71 cos 29 + 72 sin 29) - 7^. (47) 

These results can be found in Kormann et al. (1994) and Keeton & Kochanek (1998). Below 
we first discuss the simpler case without shear and then discuss some examples of the general 
case with shear (in an arbitrary direction). 

3.1 The case without shear (7 = 0) 

For this case, we have four (two) images if the source is inside (outside) a diamond-shaped 
caustic (see Fig. 1). Again, we study the sum of signed magnifications for the four image 
case. Recall that for an elliptical potential with the same power-law slope (/? = 1), the sum 
of signed magnification is exactly 2. Now for the isothermal elliptical density distribution, 
surprisingly, we find a value very close to 2.8, nearly independent of the ellipticity. For axial 
ratio in the range of 0.1 < g < 1, the sum of signed magnification varies little, from 2.65 to 
2.8. This can be demonstrated by considering the special case where the source is located 
at the origin. The lens equation can be easily solved and the total signed magnification is 
given by 

4 2 2 

Y. f^iPi = -\ 77 — ^)— + -, 77 — T^^' (48) 

^ l-u/ tan ^u I -vj tanh v 

where u = [2e/(l — e)]^/^ and v = [2e/(l + e)]^/^. The above expression can be expanded in 

series around e = 0, 

^ 14 32 2 490272 , ^, ., 

^'^^ 5 875 21896875 ^ ^ ^ ^ 

The expansion only involves even powers of e and the coefficients are very small, which 

explain why the sum is nearly independent of e (or q). Numerically, we find that the sum 

varies < 10% as long as we have four images, i.e., when the source is inside the caustics. 

To illustrate the near constancy of the sum of signed magnifications, we show in Fig. 1, 

the magnification contours for two values of q: a nearly spherical case with q = 0.9 and a 
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Figure 1. The lejt panel shows the contours for the sum of signed magnifications for q = 0.9. The axes are in units of a in 
eq. ([l4|). The thick diamond-shaped curve is the caustic. The three dashed contours (from outer to inner) have magnifications 
2.785, 2.79, and 2.795 respectively. The three solid contours (from inner to outer) have magnifications 2.8, 2.805 and 2.81, 
respectively. The outgoing contours from the cusp are those for two image configurations. The right panel is for q = 0.3. The 
dashed contour levels (from outer to inner) are 2.7, 2.725, and 2.75 respectively, while the four solid contour levels (from inner 
to outer) are 2.775, 2.8, 2.85, and 2.9, respectively. The outgoing contours (closely spaced together) from the cusp are contours 
for two image configurations. Notice the scales on the two axes are different. 

fairly elongated case with q = 0.3; these values bracket the observed light axial ratios in 
early types of galaxies. The contour shapes are qualitatively the same; the case with q = 0.9 
has a much smaller range in the total signed magnification (varies from ^ 2.79 to ~ 2.82) 
than the more elongated case. We emphasize that although the sum is nearly constant, it is 
quite different from the value (2) of the isothermal potential. This cautions against simple 
extrapolation of the magnification results obtained from elliptical potentials to elliptical 
density distributions. 



3.2 The case with shear 

The lens equation with a shear is more complicated and can be studied only numerically. 
However, we find that for most cases, the sum of signed magnification is still around 2-3. In 
Fig. 2, we show one example for 7 = 0.1, 0s = 45° and q = 0.8. For this set of parameters, 
the caustic structure is a deformed diamond, which divides the four image region (inside the 
diamond) and the two image region (outside the diamond). The sum of signed magnification 
at the origin is 2.47, while most areas inside the caustics have magnifications between 2.2 
to 2.8. 

However, there are significant exceptions to this. For some combinations of shear and 
ellipticity, the total signed magnification falls significantly outside the range of 2-3. The rea- 
son for this is that qualitatively different image configurations appear. More specifically, for 
some combinations of the 7 and e, the caustics become small, and higher-order singularities 
(such as swallowtails) emerge. Fig. 3 shows how the caustics change for q = 0.6 when the 
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Fig ure 2. Contours for the sum of signed magnifications for q = 0.9,7 = 0-1 S'lid 4'^ = 45°. The axes are in units of a in eq. 
mm. The thick diamond-shaped curve is the caustics. The three dotted contours have total magnifications 2.2, 2.3, and 2.4 
respectively. The three solid contours have magnifications 2.5, 2.6 and 2.7, respectively. The outgoing lobes originating from 
the cusp are contours for two image configurations. 

strength of shear changes but keep its angle fixed at (/)s = 90°. When the shear is <^ 0.1, 
the caustics form a diamond, similar to the case without shear. When the shear approaches 
a critical value (~ 0.15 for our case), swallowtails begin to develop, first along the y axis 
(see the top right panel), then along the x axis (lower right). For some intermediate shear 
value (7 = 0.165, lower left), complex caustics appear. Finally, when the shear becomes 
sufficiently large {^ 0.21), the caustics (not shown) again covert to a simple diamond. For 
a source at the origin, the maximum number of images is eight, as can be seen for the case 
with 7 = 0.165 (lower left panel). For this case, the total signed magnification the four 
image configurations spans a considerable range, from ?^ — 20 for various 4-image caustic 
regions. These values are considerably different from 2 found for the isothermal potential, or 
2.8 for the isothermal elliptical density distribution without shear. The magnification for the 
central region with 8 images is approximately 5.06; this magnification does seem remarkably 
constant within 1% or so. 

We have also checked for higher moments of the image positions weighted by the signed 
magnifications for the elliptical density distributions. We find that in general these moments 
no longer satisfy the simple relations we presented in §2, so we do not discuss these relations 
further. 
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Figure 3. Caustic structures for an elliptical density distribution with q = 0.6 and (pB = 90° as a function of the magnitude 
of shear 7. The axes are in units of a in eq. (|44[). The value of shear 7 is indicated in each panel. Notice how swallowtail 
singularities appear and disappear as the shear changes. For 7 > 0.2, the caustics (not shown) are again a diamond, similar to 
the one shown at the top left, only elongated along the y-axis. For 7 = 0.165, the contours for the total signed magnification 
are shown for levels of 1.5 (dotted), 3 (dashed), 5.09 (thin solid), and 8 (thick solid). The thick complex-shaped curve is the 
caustic. 

4 SUMMARY AND DISCUSSION 

In this paper, we have studied the eUiptical potentials using a new moment approach. We 
have shown that the image positions weighted by signed magnifications satisfy fairly simple 
relations for this class of potentials. For example, the total signed magnifications is 2 for an 
isothermal elliptical potential, independent of the axial ratio (ellipticity); similarly simple 
results hold for higher moments. Our results confirm and extend those of Dalai (1998) on 
the power-law elliptical potentials. These results provide insights into the lensing properties 
of elliptical potentials. With improved image positions and flux ratios for lenses, the mo- 
ment approach should allow us to understand more easily whether the power-law elliptical 
potentials can fit the observations data at all, without much numerical effort. 

However, we also find that these relations do not hold for its cousin - the singular isother- 
mal elliptical density distribution. For this elliptical density distribution without shear, the 
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total signed magnification is 2.8, again nearly independent of the ellipticity, but significantly 
different from the value (2) for an isothermal potential. For the case with shear, the mag- 
nification relations become more complicated. For most cases we checked the total signed 
magnification is still around 2-3. However, the magnification no longer changes with the 
shear quadratically as in eq. (^2]), instead it is much more sensitive. There are also sig- 
nificant exceptions for the total signed magnification for some combinations of shear and 
ellipticity. The most striking feature is that 6 or 8 image configurations can occur (see Fig. 
3). Notice, however, that the caustic structures are very small, so such image configurations 
may have too small cross-sections to be seen unless one observes a large number of quasars. 
This is likely to be true even after one takes into account the high magnification bias ex- 
pected for these cases. Such cases, if discovered, will be easy to identify since the 6 or 8 
images are located approximately on a circle. 

The comparison between magnification relations in elliptical density distribution and 
elliptical potentials highlights the fact that magnifications are highly sensitive to the detailed 
form one uses for the lensing potential, much more so than the image positions; this is 
easily understood intuitively since magnifications involve higher order derivatives of the 
lensing potential which tend to amplify small differences in the parameters. Since even the 
more realistic elliptical density distribution is only an approximation for the real lensing 
potential, as evidenced by the twisted isophotes and many other irregular features in the 
light distribution, it is a reasonable extrapolation that the predicted magnifications from 
simple lens models can be significantly different from the observed ones. Indeed, for some 
radio lenses, one can show that even substructures with mass ^ IO^Mq can significantly 
influence the magnifications (Mao & Schneider 1998). It means that it will be difficult 
to model the flux ratios using simple lens models. However, this also means that with the 
steadily improving lens flux ratios and astrometries (e.g., see Blanton, Turner & Wambsganss 
1998 for 2237+0305), we may be able to probe the structures in lensing galaxies in great 
detail and discriminate different classes of lensing potentials with confidence. 
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